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We consider the interaction of an harmonic oscillator with the quantum field via radiation 
pressure. We show that a 'Schrodinger cat' state decoheres in a time scale that depends on the degree 
of 'classicality' of the state components, and which may be much shorter than the relaxation time 
scale associated to the dynamical Casimir effect. We also show that decoherence is a consequence of 
the entanglement between the quantum states of the oscillator and field two-photon states. With the 
help of the fluctuation-dissipation theorem, we derive a relation between decoherence and damping 

' rates valid for arbitrary values of the temperature of the field. Coherent states are selected by the 

' interaction as pointer states. 

o : 

(N ■ 

'. I. INTRODUCTION 

• Superposition states have an important role in the formalism of quantum mechanics. However, they are in flagrant 
contradiction with our classical world when the components correspond to macroscopically distinguishable states. 
T-H , The reason why these states are not encountered round the corner is decoherence, a process by which the interaction 
between the degrees of freedom of the system in question with any other degrees of freedom, either internal or external 
(the so-called environment), leads to a suppression of the coherence between the components of the superposition [Q. 
Even if this coupling is very weak, the decoherence rate may be huge, resulting in a very fast decay of these "weird" 
states and in the emergence of the classical world. Recent developments in technology now allow one to study in 
r real-time the process of decoherence in the lab. For example, over the past several years techniques have been 

■ developed to generate mesoscopic superpositions of motional states of trapped ions 1^ , and of photon states in cavity 
, quantum electrodynamics In these cases decoherence due to the coupling with the ambient reservoirs was observed, 

■ confirming the expectation that the decoherence rate is faster, the larger and more separated the state components 
, are Q|. Recently another experiment has succeeded in "engineering" the environment in the context of trapped ions, 

^ ■ studying scaling laws of decoherence theory for a variety of reservoirs in a wide range of parameters . 
CIh! Usually, decoherence is analyzed in the framework of heuristic models that describe phenomenological dissipation 
^l) i (the reservoir is taken to be a collection of harmonic oscillators, coupled linearly to the position operator of the 
^ ' system [^0). In this paper, we consider instead an ab initio model for decoherence of a particle in a harmonic 
potential, scattering the radiation field (at temperature T), which then plays the role of the reservoir. Starting 
from first principles, we show that the resulting radiation pressure coupling with the field transforms an initial pure 
superposition state of the particle into a statistical mixture. 

Of special relevance is the limit T = 0. In this case, the reservoir is the quantum vacuum field, which dissipates the 
mechanical energy of the oscillating particle (or 'mirror'). This effect is associated to the emission of pairs of photons, 
the so-called dynamical Casimir effect. Much work has been done on quantum radiation from moving mirrors [||. 
i Important properties like the spectrum of the emitted radiation the time evolution of the energy- momentum 
tensor [ p^ , the total radiated intensity and the dissipative radiation pressure on the particle (radiation reaction force 
corresponding to the photon emission effect) ||ll| have been considered. Here we focus our attention on the 
particle as the system of interest, and show that decoherence is a consequence of the entanglement between particle 
and field two-photon states. This result has fundamental implications, for it shows that any particle not transparent 
to the radiation is unavoidably under the action of decoherence through the radiation pressure coupling with vacuum 
fluctuations. 

The zero temperature limit was briefly discussed in our previous letter jisf . This article presents results for finite 
values of temperature, as well as a detailed discussion of the case T = 0. The formalism relies on the ID scalar model 
for the field, but extensions to 3D electromagnetic field are also discussed, allowing us to address the question of orders 
of magnitude. The paper is organized as follows. In section II we start from the Hamiltonian model for the radiation 
pressure coupling, and then derive a master equation for the particle. In section III we discuss how the environment 
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selects a prefered basis in the particle's Hilbert space, the pointer basis. In Section IV we derive a general relation 
between decoherence and damping rates at arbitrary temperature by means of the fluctuation-dissipation theorem. 
The zero and high temperature limits, including extensions to the 3D electromagnetic model, are discussed in Sections 
V and VI. Section VII contains our conclusions. Finally, in the appendix an alternative, simpler derivation of the 
decoherence rate is given, which is based on the entanglement between the particle and two-photon states. 



II. THE MODEL 



Most treatments of the dynamical Casimir effect are based on the assumption that the mirror follows a prescribed 
trajectory, thus neglecting the recoil effect. However, in this paper we want to focus on the mirror as a dynamical 
quantum system, hence the need to tackle the full mirror-plus-field dynamics. This has already been addressed in 
the framework of linear response theory in order to calculate the fluctuations of the position of a dispersive mirror 
driven by vacuum radiation pressure |Q, and related calculations have been performed in Refs. p^-p^ to derive 
mass corrections caused by the interaction with the field. 

We consider a nonrelativistic partially reflecting mirror of mass M (with position q and momentum p) in a harmonic 
potential of frequency cjo, and under the action of radiation pressure. We take a scalar field in 1-1-1 dimensions, which 
mimics the electromagnetic field modes that propagate along the direction perpendicular to the plane of the mirror. 
Extensions to the real 3-1-1 case are analyzed in Sees. V and VI. We neglect third and higher order terms in w/c, 
where v is the mirror's velocity (we set c — \ hereafter, except when an explicit evaluation of orders of magnitude is 



required). Our starting point is the Hamiltonian formalism developed in Refs. |16| and |17| (although these references 



consider a free mirror, the extension to the harmonic oscillator is straightforward). The total Hamiltonian is 

H = Hm + Hp + Hint, (1) 



where 



is the harmonic oscillator Hamiltonian for the mirror, and 

HF = J^[n' + {d.,cf>r] + n^\x = 0) (3) 

is the free Hamiltonian for the field (j) and its momentum canonically conjugated H = dt4>. The second term in the 
r.-h.-s. of Eq. (^) is associated to the boundary condition of a partially-reflecting mirror at rest at a; = 0. In the 
context of the plasma sheet model of Ref. [|l^, it corresponds to the kinetic energy of the plasma charged particles. 
The coupling constant plays the role of a transparency frequency, since from Eq. (^) one derives the boundary 
condition 

5,</.(0+) - 5,(/.(0-) = 2r!0(O) 
{(j) is continuous at x = 0), which yields a frequency-dependent reflection amplitude ||lj] |p^ : 

R{uj) = -I (4) 

Finally, the interaction Hamiltonian is given by 

Hint = -— + — - -Q-ch^ix = 0)^, (5) 

where V — — J dxdxff) dt4' is the field momentum operator. Hint describes, to second order in v/c, the modification 
of the boundary condition for the field due to the motion of the mirror, which in its turn is affected by the field 
radiation pressure. Thus, it provides a coupling between the harmonic oscillator and the field, to be treated within 
perturbation theory. The small perturbation parameter is v/c, and not the transparency frequency fl, which may be 
arbitrarily large. The first term in Eq. (|^) is responsible for the effect of decoherence to be discussed here. It also 
accounts for the effects of emission of photons, dissipation of the mirror's energy, and part of the mass correction. 

We calculate the density matrix i5{t) of the combined mirror-plus-field system using second order perturbation 
theory, and trace over the field operators to derive the master equation for the mirror's density matrix p(t) |l^. We 
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assume that at t = the mirror and field are not correlated: /5(0) = p(0) (S) Pf, where is the density matrix of the 
field (assumed to be in some steady state; later in this section we take a thermal equilibrium state). We find 

tnm = [Hm, p{t)] ^^-^^ \p\ pm (6) 



2hM 







where the superscript / indicates the operators to be taken in the interaction picture. The second term in the r.-h.-s. 
of Eq. ^ is the contribution in first-order of perturbation theory of the term in the interaction Hamiltonian (see 
Eq. (^). It corresponds to a (cut-off dependent) mass correction given by 

AMi = 17(02(0)), (7) 

as already found in Refs. [ p^ and |^6|. The (anti-)symmetric second order correlation function (^) a is defined as 

a{t) = C{t) + C{-t), (8) 

m = C{t) - Ci-t), (9) 

with 

C{t)^{V'{t)V'{Q))-{Vf. (10) 

When computing the correlation functions, we take the unperturbed field, which corresponds to the static boundary 
condition (eigenfunctions of Hp). 

Replacing the free evolution for p^ {—t') in (^) yields 

ihp = [Hm - -jf^j^^p] - r(t)b,{g,p}] - ^^i(Ob, [p,p]] ~ ^^2(t)[p, hp]]- (n) 

The total mass correction is AM — AAfi + AM2, where AM2, as well as the remaining coefficients in (pT|), originate 
from the first term in the r.-h.-s. of Eq. (^, taken in second-order perturbation theory. Their meanings are best 
understood when writing the Fokker-Planck equation for the Wigner function W{x,p, t) : 



dtW = -(1 - AM/M)4-:d^W + MLolxd^W + 2VdJxW) + Di — W - D2Tr^W. (12) 



1 

AM2 and the damping coefficient T are calculated from the anti-symmetric correlation function 



AKhit)^- dt'cos{ujot')^{t'), (13) 
Jo 
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whereas the diffusion coefficients are associated to the symmetric correlation function: 
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Di{t)^^ J^dt'co8{u;ot')<j(t'), (15) 

D2{t) = ^J^ dt' sin(c.o<') <jit'). (16) 

We assume the field to be in a thermal state (temperature T), and take the following strategy to calculate the 
momentum correlation functions. The time derivative of the field momentum is minus the radiation pressure force on 
the mirror iflql: 
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— -2r!0(o,i)a,0(o,i), (17) 

where dx(f>{0,t) = [dx(f>{0^) + dx(f>{0~)]/2. Using Eq. (0), we calculate C{t) by integrating the correlation function of 
the field calculated at x = : 

t 



c{t) = {2nf dh [(0(o,ii)40(o,ii)(/)(o,t2)4</'(o,i2))- 

J —OO J —OQ 

(0(0, ti)dxcb{0, ii))(0(O, i2)40(O, t2))] (18) 



The equal-time second-order correlation function in ( |18| ) corresponds to the force on the static (single) mirror. It 
vanishes since the vacuum radiation pressures exerted on each side of the mirror are in perfect equilibrium. On the 
other hand, the fourth-order correlation function may be expressed as a sum of second order correlation functions 
(with the fields taken at different times), which are calculated with the help of the normal mode expansion for the 
field. They are directly connected to the average number of photons riuj = l/[exp(?ia;/r) — 1] in the mode of frequency 
uj at temperature T (we take the Boltzmann constant fee = 1). It is useful to write the result in the Fourier domain, 
the Fourier transform of the anti-symmetric correlation function ^{t) being defined as 



CM = J dtexp{iLut)^{t). (19) 

Eqs. (|) and (|l|) yield 

^[u;]=e[c^]+eH (20) 

where 

= {2/Tr)h^nciuj/n), (21) 

with 

C(it) = ln(l + ■ti^)/(2-u) + (arctanu)/!*^ - 1/u, (22) 
represents the correlation function at T = (vacuum fluctuations), whereas 



with 

G{lo,lo') = \uj' - uo\ {n\^,_^\ - €{lo' - Lo)n^,) , (24) 

represents the thermal fluctuations (e is the sign function). 

Symmetric and anti-symmetric correlation functions for a system in thermal equilibrium are related in a very general 
way H pi: 



tanh( 



2T' 



According to Eqs. ( |14| ) and (^_5|), this result provides a general relation between diffusion and damping, in the spirit of 
the fluctuation-dissipation theorem. This relation is particularly simple for the asymptotic values of the coefficients 
r(t) and Di{t) at i — > oo. Since the integrands in Eqs. ( p^ and (^5|) are even functions of time, we may extend the 
integration range to — oo, yielding 

- CN] (26) 
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and 



4 



Thus, the asymptotic values of T and Di are directly connected to the fluctuations at the mechanical frequency luq, 
allowing us to derive, from Eq. (|2^), a simple and general relation between these two coefficients. On the other hand, 
no such simple connection exists for the remaining time-dependent coefficients, AM2 and -D2, whose asymptotic values 
result from the joint contribution of the whole spectrum of fluctuations fl^ . 
Combining Eqs. (H)-®, we find 

Di - — , (28) 

MtJotanh(^) ^ ' 



a clear manifestation of the fluctuation-dissipation theorem. According to Eq. (|2g), the temperature dependence of 
the diffusion coefficient is determined, apart from the T dependence of the damping coefficient F (to be discussed 
later), by the relative importance of the thermal fluctuations (and their corresponding energy T) with respect to 
quantum fluctuations (and their corresponding zero-point energy ?iujq/2). In the high-temperature limit, T ^ ?ia;o/2, 
Eq. ( p^ ) yields Di = 2Tr /{Mluq). In the theory of Brownian motion, F is a T independent phenomenological constant, 
and hence the diffusion coefficient is a linear function of temperature in this limit. Here, however, F has an explicit 
temperature dependence, to be analyzed in Sec. V. 

From Eq. (|28|), we shall derive a relation between decoherence and damping time scales, valid for any temperature 
T. Before considering a specific superposition state, however, we discuss, in the next section, the degree of sensitivity 
of different states in the Hilbert space to the action of decoherence. We also analyze in more detail the precise meaning 
of t —^ 00 (in the par ticular case of T = 0), in order to know how fast the coefficients approach their asymptotic 
values. From Eqs. (|lj) and (|l^) alone it may be shown, in a general way, that a sufficient condition is t ^ l/^o, but 
in some cases the convergence may be much faster. 

III. POINTER STATES 

Different criteria have been introduced in the literature in order to find out the states in the Hilbert space that are 
most robust under the interaction with the environment and behave more classically [p"9|-^ . Here we shall follow the 
one introduced by Zurek, the so-called "predictability sieve" . The idea is to take every possible state of the Hilbert 
space, calculate its entropy, and order the states in a tower according to increasing entropy. The most classical states 
are those that lie at the bottom of that tower, and correspond to the most predictable ones. For these 'pointer' states, 
information loss due to the interaction with the environment is minimal. This philosophy is put in quantitative terms 
by minimizing the linear entropy of the system, 

s[p] = l- trp2, (29) 

which is zero for a pure state and greater than zero for a statistical mixture. In general this is a difficult problem 
because complicated entanglement between system and environment develops on account of their mutual interaction. 
So far, results have been successfully derived assuming that the initial state of the system is pure. Here we follow the 
same approach, and calculate the rate of entropy increase starting from the master equation ([Tl|). We assume that 
the state is nearly pure at time t to find 

m = - 1) + ^Ap^ + (30) 

where Ap^ = (jp) — (p)'^ is the momentum dispersion and (Tg_p = p}) — 2{p){q). Here (...) — tr(. ..p), and 
all operators are evaluated at the same time t. The first term in ( pO|) leads to a decrease of entropy s{t) = 1 — 
exp(2 J^r{t')dt'), he nee damping tries to localize the state competing against diffusion. This decrease is independent 
of the initial state of the system, and therefore is irrelevant for determining the pointer states. 

We assume that the typical decoherence time scale is much larger than the period of free oscillation 2tt/ujo, so that 
we may integrate Eq. (|30| ) to find the entropy at an intermediate time t = n2TT/u!o. We take n 3> 1, allowing us 
to replace the time dependent coefficients by their constant asymptotic values, but assume that t is much shorter 
than the decoherence time scale, in order to be consistent with the small-entropy approximation underlying Eq. ([3C|). 
Moreover, in this weak coupling limit, we may take the free evolution (corresponding to the harmonic oscillator 
Hamiltonian Hm) for the mirror's operators in Eq. (|30|). The correlation function cr^.p oscillates around zero, and 
then does not contribute to s(t), whereas the free evolution of Ap^(t) mixes up position and momentum fluctuations, 
yielding 

s(t) = 2t^ [iAp)l + {Muj^f{Aq)l ~ Mhuj^] , (31) 
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where (Ap)q and {Aq)^ represent the dispersions for the initial state. From Eq. (pT|), we find that the minimum 
entropy given the constraint AqAp > h/2 is for Aq^ = h/ {2Mujo), Ap^ = MTiloq/2. Hence, as in the Caldeira-Leggett 
model, and for any temperature of the field, the pointer basis consists of coherent states [|o|. In this weak coupling 
approximation, the minimum value corresponds to sir) — 0, hence the increase of entropy of a coherent state is a 
higher order effect. 
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D,(t)/D, 

D,(t) / (M a)„ Dj) 

r(t) / r 



FIG. 1. Diffusion and damping coefficients for zero temperature as functions of time in the perfectly-reflecting limit, 
LJo/i^ = 10"* < 1. Here Di = h'^iJo/12nM^ and T = ?iaJo/127rA'/ are the asymptotic limits of Di{t) and r(t). The in- 
sets show the behavior of these two time-dependent coefficients for short times. 

The crucial approximation in the derivation of (|^) from ( |30| ) is the replacement of the time dependent coefficients by 
their finite, constant asymptotic values. It is instructive to analyze in more detail the behavior of the coefficients and 
its connection with entropy production. As an example, we take T = 0, and consider first the perfectly-refiecting limit, 
which corresponds to ^ i^, for in this case the relevant field modes have frequencies much smaller than the mirror's 
transparency frequency. In Fig. 1 we plot the diffusion and damping coefficients as functions of ujot for ajo/n = 10~^ 
and T = 0. The damping coefficient F approaches its asymptotic value very fast, for t ~ whereas Di{t) develops 
an initial jolt for times of the order of fi"^ and then decreases to the asymptotic value (I?i)perf = ft^wo/(127rM^) for 
t ^ l/ujQ. When we integrate Eq. ( ^0| ) over many periods of oscillation, the contribution to the entropy of the initial 
jolt is negligible, allowing us to replace Di by its asymptotic value. 



D,(t)/D, 

D,(t)/(Mco„D,) 

r(t)/r 




FIG. 2. Diffusion and damping coefficients for zero temperature T = as a function of time in the high-transmission limit 
Loo/^ = 10* > 1. Here Di = h'^Q.'^ ln(ajo/n)/27rAf^a;o and T = TiQ'^ ln{Loo/^) /2-n-M are the asymptotic limffs of Di{t) and r{t). 
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In Ref. IgTj, it was shown that no net entropy is produced for the Caldeira-Leggett model with an adiabatic 
environment, since aU the time-dependent coefhcients are oscillating functions around a zero mean. At first sight, the 
same would happen in our model when ojq 3> ^, for in this case the dominant field frequencies are slow with respect 
to the mirror's translational time scale. However, as discussed in Section V, the spectral density decays too 

slowly for lu Q, and as a consequence field frequencies of the order of ujq provide a significant contribution even 
in this limit. Thus, one cannot ascribe a frequency cut-off to the environment such that the typical frequency of the 
system uq is much greater than the maximum frequency of the environment. Therefore, the vacuum field does not 
behave adiabatically in the sense of In our case instead, the diffusion coefficients oscillate around a non-zero 

value, leading to a net entropy increase. In Fig. 2, we plot the diffusion and damping coefficients as functions of ujot 
for wo/rj = 10^ and T = 0. They oscillate around their asymptotic values with (angular) frequency ojq and with an 
amplitude of oscillation that decays in a time t ~ 1/Q ^ I/wq. The oscillatory terms do not contribute to the entropy 
increase when we average over many oscillations. Hence Eq. (^) also holds in this case, although the rate of entropy 
increase is much smaller than in the perfectly-reflecting limit. 



IV. DECOHERENCE VERSUS DAMPING 



In this section, we derive a general relation between damping and decoherence time scales, starting from the 
fluctuation-dissipation result given by Eq. (|2^). As an extreme case of decoherent dynamics, we consider a superposi- 
tion of two coherent states, since they correspond to the pointer states, according to the results of Sec. III. Specifically, 
we take at i = the even superposition state |V')o = + | — Oi.))/^/2, with a — iPo/^/2MhIjo, so that the coherent 
states are initially along the momentum axis in phase space, and ±Pq are the average values of momentum of the 
components at t = 0. We also assume that \a\ 3> 1, hence the average energy of the state components is much larger 
than the zero-point energy. The corresponding Wigner function is 



1 

W ^ Wrn + —r exp 

ttTi 



7 7 



2(Ago)2 2(Apo) 



cos(^), (32) 



where Ago = yhJ{2MujQ) and Apo = ?i/(2Ago) are the position and momentum uncertainties of the ground state. 
Wm corresponds to the statistical mixture 

p,n^[\a){a\ + \-a){-a\)/2. (33) 

In phase space, Wm has two Gaussian peaks along the momentum axis at ±Po- Pm is a classical state in the sense that 
Wrm being positive defined, may be interpreted as a probability distribution in phase space. On the other hand, the 
nonclassical nature of the superposition state is featured by the remaining term in Eq. ( |34 ) , representing the coherent 
interference between the two state components, and which oscillates into negative values along the position axis. 

Diffusion along position, associated to the coefficient Di , averages out the oscillations of the interference term at a 
rate to be calculated from the Fokker-Planck equation ( p^ | |2^ . According to Eq. (32), the oscillations are faster 



the higher the value of Pq, so that td is a decreasing function of \a\. As in Section HI, we assume that decoherence 
is very slow, 1/td «C wq, so that several free oscillations take place before coherence is lost. In this limit, the particle 
has enough time to probe the harmonic potential before diffusion takes place, and as a consequence decoherence is 
governed by the asymptotic value of which is directly connected to the field fiucutuations at the frequency of 
oscillation wo, according to Eq. (|2^). This condition holds for most experiments, where mesoscopic superpositions 
are employed so as to render decoherence slow enough to be measured Q ||]. Moreover, it always applies in the case 
of vacuum radiation pressure (T = 0), as shown in Sec. V. Diffusion is maximum when the state components are 
along the momentum axis: from (|3^), we find d^W ~ —{2PQ/Kf'W; and vanishes when the two wavepackets reach 
the turning points in the harmonic potential. The average over many oscillation yields 

1 _ 1 (dlw\ _2PlD, 

\ / max 

that combined with the fluctuation-dissipation theorem, Eq. (^), yields the following result for the decoherence time 
td : 

^<i = ^t^-h(^)^. (35) 



7 



AT = (or more generally, for T <C fiuJo), Eq. (35) yields = l/(4|apr). This result may be written in terms of the 
distance AP = 2Po between the two components in phase space at t = 0, or in terms of the distance AQ = AP/ (Mluq) 
at t — t:/2luo : 



The interpretation of ( |3^ ) is clear: decoherence is faster, the more separated the state components in phase space 
are. Here the zero-point fluctuations define the reference of distance in phase space. At high temperatures, on the 
other hand, this reference is provided by the thermal de Broglie wavelength At — Ti/\/2MT. In fact, ( ^ ) yields, for 

. _hcJo 1 1 [XtVi 

Eq. ( |37|) also shows that the ratio between decoherence and damping rates is larger at high temperatures by the factor 
2T/{hiJo). 

When written in terms of distances in phase space, the results above are also valid for more general superposition 
states, like (|0) + |a))/\/2. Mor cover, their range of validity is not limited to the radiation pressure coupling considered 
here. In fact, Eqs. ( p6| ) and ( ^7| ) are in perfect agreement with the results obtained in the framework of the Caldeira- 
Legget phenomenological model for quantum dissipation |^. This is hardly surprising, since they rely on general 
properties of the correlation functions associated to the fluctuation-dissipation theorem. Eq. (^5|), which interpolates 
the low and high temperature limits, is also discussed in Ref. ||^, in the context of a two- level system. The temperature 
dependence for the ratio between decoherence and damping times has a simple interpretation: at T > 0, the time 
scale for the relaxation of the populations is shorter than l/F exactly by the factor tanh(?itJo/2T), on account of 
the contribution of absorption and stimulated emission. Here this factor originates from the general relation between 
symmetric and anti-symmetric correlation functions, Eq. ([2^), which is at the heart of the fluctuation-dissipation 
theorem. 

The peculiarities of the radiation pressure model considered here are contained in the damping rate F. Rather than 
a phenomenological input parameter, it is computed from first principles, first for T = in Sec. V, and then for 
T > hcja in Sec. VI. 



V. VACUUM FIELD 

At T = 0, the spectral density is given by Eqs. (|2l|) and (p^. This result is more easily obtained from the following 
argument (a similar method, applied for the force correlation function, may be found in Refs. ||2^ and ||2^). Since 
V is quadratic in the field operators, the correlation function C{t) may be calculated from the two-photon matrix 
elements of the momentum operator as follows: 

1 I'OO I'OC 

C{t) = - duJi diU2{0\V{t)\uJi,UJ2){uJl,iO2\V\0). (38) 



2 Jo 

We have {0\V{t)\uJi,uj2) = cxp[—i{LOi^LJ2)t]{0\V{0)\iJi^uj2) since only the annihilation operators contribute, and hence 

rOO roc 

CM^Tt/ duJl duJ2\{0\Vi0)\LOl,UJ2)fSiuJ-UJl-LO2). (39) 

Jo Jo 

Thus, at r = the fluctuations at the (positive) frequency lo originate from two-photon states \u}i,uj2) such that 
LU1+UJ2 = uj- In the dynamical Casimir effect, the oscillation at the mechanical frequency loq gives rise to the emission 
of pairs of photons of frequencies loi and uj2, such that = uji + lo2- On the other hand, according to (p6[), F originates 
from the fluctuations at frequency loq, and hence 

rco POO 

r=-^y^ dWij^ duJ2\{uJu^2\VmH{uJo-UJ^~L02), (40) 

rendering explicit the connection between damping and the photon emission effect. In the appendix, we present 
an alternative derivation of (pO|), starting from the two-photon emission probabilities and making use of energy 
conservation. 
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Eq. (39) also shows that C[uj] vanishes for negative frequencies, and as a consequence, a\uj] = e(a;)^[a;] in agreement 
with Eq. (H). Finally, the result of Eq. (|l|) follows from (H) by using again Eq. (0) [Q. In Fig. 3 , we plot 
as a function of its argument. According to Eq. (|2l|), the transparency frequency fl defines a scale for the behavior 
of the spectrum of fluctuations in vacuum. For lu <ti Q, the spectrum is linear: (^{uj/fl) « a;/(6ri), and goes to zero 
slowly, as \ii{u! / fl) / (lu / Q) , for lu ^ Q, due to the high-frequency transparency of the mirror. 
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FIG. 3. Spectral density for zero temperature. 



The damping coefHcient at zero temperature is obtained from Eqs. (|2l|) and (^6|), or alternatively from Eq 



In 



the perfectly-reflecting hmit, lvq ^ 17, Eq. ( p| ) yields 



r = 



UttM 



Wo- 



(41) 



(42) 



Thus, the damping induced by the Casimir effect is a small perturbation of the free harmonic oscillations. The ratio 
between the zero-point energy and the rest energy appearing in (^2|) is also of the order of the recoil velocity of the 
mirror divided by c, which is, as explained in Sec. II, the small parameter of the perturbation approach leading to 
the master equation (|l^). For larger values of wo/fl, the damping as given by (|4^) is still smaller, since vacuum 
frequencies of the order of uq are not well reflected by the mirror in this case. 

Eq. ( |42[ ) is directly connected to the well-known formula for the dissipative Casimir force on a single perfect moving 
mirror |11| , F = hx / (Gvr) , for the equation of motion then reads pq] 



LUqX + 



hx 

GttM' 



(43) 



whose solution when hu!o/M ^ 1 is 



xoe"*""* exp 



hujQ 
127rA/' 



The decoherence time scale at T = in the perfectly-reflecting limit is derived from Eqs. (p6|) and (E^): 



3 27r 

V' UJQ 



(44) 



where v = Pq/M is the initial velocity of the wavepackets. Being of the order of (w/c)^, the decoherence rate is 
very small at T = (or, at any rate, in the nonrelativistic limit considered here). Since wo^d 3> 1, decoherence is 
the cumulative effect of several free oscillations in the harmonic well, which justifies the approach employed in the 
derivation of (^) and the use of the asymptotic value for Di{t). 

In order to further understand how the dynamical Casimir effect engenders decoherence, we present, in the ap- 
pendix, an alternative approach, where we follow the evolution of the complete oscillator-plus-field quantum state. It 
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shows that the superposition state decoheres because the two wavepacket components oscillating out-of-phase yield 
amplitudes for emission of photon pairs with opposite signs. As a consequence, an entangled mirror-plus-field state is 
developed, given by 



oc 



|«')At = S(Ai)|V')e® |0) + i / duji I dw2&K,'^2;Ai) |V')o«) ki,t^2), (45) 

^ Jo Jo 

where \il')o = — | — ct))/V2 is the odd superposition state, 0^2; Ai) is the amplitude for emission of a photon 
pair with frequencies loi and uj2 during At (the explicit expressions are given in the appendix), and B is the amplitude 
for persistence in the vacuum state: 

|B(Ai)|2 = l-- / dcol / dcj2\b{cji,oj2,m^- (46) 
^ Jo Jo 

The density operators of the odd and even superposition states differ by the sign of the interference term, pint — 
P — Pm {Pm is defined in Eq. (|33|)). Accordingly, when computing the reduced density matrix of the mirror, p{At) = 
trficid(|^')At(2lD: the contribution of the two-photon states in Eq. ( p5[ ) reduces the coherence of the state. With the 
help of Eq. (Eq) , we find 

Apint = Pint (Ai) -Pint (0) = -- / duJi | &(^1 , ^-'2 ; Ai) | ^ pint (0) . (47) 

^ Jo Jo 

The two-photon probabilities are proportional to At and connected to the damping rate T as discussed in the appendix. 
Hence Eq. (|^ yields 

dpint _ Apint Pint 

-dt ^ - "17' ^ ^ 

with id given by Eq. (|3^). 

In this derivation, the expression for the emission amplitudes per se are not of any relevance — only its connection 
with the damping rate T is important. This connection is based on the principle of energy conservation: the energy 
of the oscillator is damped at the rate at which energy is radiated. Since this argument also holds for the real 
3-1-1 electromagnetic field, we may extend our results by replacing the 3D result for F into Eq. (^. The dissipative 
dynamical Casimir force on an oscillating (frequency c^o) perfectly-refiecting sphere was obtained in Rcf. [ p3[ . Usually, 
the sphere is very small when compared with the wavelength of the relevant vacuum fluctuations, which are of the 
order of 27r/uJo. When ojqR ^ 1, the force on the sphere of radius R is given by 

_ (9) /4g^ 

" 6487r ' 

where a;*^^^ is the ninth time derivative of the position of the sphere. Following again the method of Eq. (^3|), we 
calculate the damping rate F from the equation of motion. We find 

^ - 1296^ M ^^^> 

showing that the coupling with the vacuum field is reduced, as compared with the ID case, by the (very small) factor 
{ujqR)^. This reduction factor accounts for the inefficient coupling of the small particle, which scatters field modes of 
very long wavelengths. Using Eq. (^6|), we find that the decoherence time increases by the same factor: 

324 1 27r 

*d = — 7 — ^ — , 51) 

hence decoherence through radiation pressure is a tiny effect at T = 0. At finite temperatures, however, the effect 
may be significant, as discussed in the next section. 

VI. HIGH-TEMPERATURE LIMIT 

In this section, we compute the damping and decoherence rates when T ^ Hloo- In this limit, vacuum fluctuations 
are negligible when compared with thermal fluctuations, and the dominant contribution in Eq. (20) comes from 



10 



which is given by Eq. (^3|). When the temperature is also higher than the cut-ofF energy hi}, all relevant frequencies 
in (p3|), which are smaller or of the order of ft, are much smaller than T/h. Then, we may take the approximation 
riui' « T/{hco'), yielding 

e[uo]^2^. (52) 

UJQ 



Replacing (|52D into (|26|) yields 



in agreement with the result for the viscous radiation pressure force obtained in Ref. F = —flTq{t). 

From a practical point-of-view, the opposite limit, hujo <^ T <^ hfl is more interesting for particles that scatter 
visible light {il ~ lO^^Hz). In this case, the corresponding reflectivity amplitude is approximately constant for 

the field modes whose frequencies are smaller or of the order of T/h. As a consequence, we may neglect the Lorentzian 
fall-off in (|23|). Moreover, we replace the thermal photon number rii^'-i^g in (pj) by 



n^'-^,~[exp{hu;'/T){l-hiUo/T)-l] \ (54) 
Neglecting second and higher order terms in hujo/T, wc find 

G(^o,w) = — 7 — (55) 



From (EST) and (55) we find 



e[u:a] = ^ — , (56) 



and then 



t: T' 



2 



which is also in agreement with Ref. [Q. It corresponds to the high-temperature, perfectly-reflecting limit. Here T 
plays the role of frequency cut-off instead of f2, so that the resulting damping rate is independent of the latter. 

The dissipative force in the high temperature limit may be interpreted as the effect of Doppler shift of the reflected 
thermal photons |2^. For a photon of frequency w, the frequency shift is Aw = ±.2ujq, the plus and minus signs 
applying for counter and co-propagating cases. Hence the motion gives rise to an unbalance between the radiation 
pressure exerted on each side of the mirror, corresponding to a momentum transfer AP = 2 Ai? q, where AiJ is the 
reflected energy during a time interval At. In terms of the density of modes g(w), we have 

/•oo 

AE= duj\R{uj)\'^g{uj)n^hLu, (58) 

"'0 

where |i?(ci;)p represents the mirror reflectivity (the reflection amplitude R is given by Eq. (|4[)). From AE, the friction 
force is obtained through 

AE , , 

F - -2—q. (59) 

In the ID case, the density of modes is frequency independent: g{Lo)duj = {L/Tr)doj, where L = At is the length 
of the quantization box. When replaced into Eq. (|5^), this result leads, with the help of (|59| ), to expressions for the 
force in agreement with our results for F, except for a factor of 2 when hujQ ^ T <C hil, [E9L In the 3D case, on the 
other hand, we have 

V 

g(uj)duj = —oj dio, (60) 

where V = AAt is the quantization volume, A being the surface of the mirror (in this case, for simplicity, we assume 
a flat rather than spherical mirror). In the limit huji^ <C T ^ hVl, Eqs. (psh and (|60|) yield 
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Af^X exp(ac^/r)-l " 15^- ^ ' 

As expected, the reflected power features the T'^ dependence of the Stefan-Bohzmann law, since it is proportional to 
the total blackbody radiation energy in this limit. The friction force is found by replacing Eq. ( |6l| ) into (59). The 
resulting damping coefRcient is given by 

Since we have neglected diffraction at the borders of the mirror, this result only applies when the mirror is much 
larger than the thermal photon wavelength Ath = IttTi/T. 

The decoherence time is then found by replacing ( p^ into (p7|), which connects damping and decoherence in the 
high-temperature limit (we re-introduce the speed of light c in order to allow an evaluation of the orders of magnitude): 

As a numerical example, we take T = bQK, which gives Ath = 2.9 X lO-^'m, and A = 1 mm . In this case, diffraction 
is negligible, and Eq. ( |6^ ) yields id[s] = 1.0 x 10~^''/(A<5^[ni^]), showing that decoherence is very fast even when the 
distance between the wavepackets is, for instance, in the nanometer range — in this case the decoherence time is of 
the order of a micro-second. Since td scales as it is still shorter, by a factor ^ 8 x 10'^, at room temperature. 

Note, however, that Eq. (^3|) only applies when wo^d ^ 1, the basic assumption that allowed us to relate decoherence 
and damping time scales with the help of the fluctuation-dissipation theorem pOl . 



VII. CONCLUSIONS 



As in the Caldeira-Legget model [g^ , coherent states are the most robust when the radiation pressure coupling with 
the quantum field is considered. This is amazingly in line with their well-known status of 'quasi-classical' states, i.e., 
the closest possible representation of a classical oscillation in a harmonic potential well. Superpositions of coherent 
states decay into a mixture at a rate proportional to the damping rate and to the squared distance in phase space. 
The ratio between decoherence and damping rates is a simple hyperbolic increasing function of temperature, which 
interpolates the zero and high temperature limits. It originates from the general relation between symmetric and 
anti-symmetric correlation functions, associated to the fluctuation-dissipation theorem. Thus, the particular nature 
of the model for the coupling with the reservoir seems to be immaterial, as far as the connection between damping and 
decoherence is concerned. Note, however, that the validity of this result is limited by the assumption that decoherence 
is slow compared to the free oscillations. 

We have shown that the radiation pressure exerted by thermal photons is a very efficient source of decoherence, 
although the corresponding energy damping effect, associated to the Doppler frequency shift of the reflected photons, 
is usually negligible. At T = 0, the energy damping is associated to the emission of photon pairs (dynamical Casimir 
effect). The dominant contribution comes from vacuum fluctuations corresponding to wavelengths of the order of 
27rc/a'o, which is usually much greater than the size of the oscillator. As a consequence, the radiation pressure 
coupling is inefhcient, and both damping and decoherence rates become very small. It is however remarkable, from 
a theoretical point-of-view, that the mere inclusion of an unavoidable, intrinsically quantum effect, is sufficient (in 
principle) to engender decoherence, and by that means restoring, although in a very long time scale, the classical 
world. 

We are grateful to A. Calogeracos and G. Barton for correspondence, and to J. Dziarmaga, A. Lambrecht, M.-T. 
Jaekel and S. Reynaud for discussions. P. A. M. N. thanks CNPq, PRONEX and FAPERJ for partial financial 
support. 



APPENDIX A: ENTANGLEMENT WITH TWO-PHOTON STATES 



In this appendix, we present an alternative, simpler derivation of the decoherence time scale at T = 0, which shows 
more clearly how the dynamical Casimir effect modifies the quantum phase of a superposition state and engenders 
decoherence. Instead of tracing over the field, we follow its evolution during many periods of free oscillation. We first 
take, at t = 0, the mirror-plus-field state |a) (E) |0) (|0) represents the vacuum field state), where \a) is a coherent 
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state of large amplitude: \a\ ^ 1. We take a = i q{0) y/ M / 2hujo , so that \a) is a 'semiclassical' state associated to 
a minimum uncertainty wave-packet whose initial velocity is q{0). We have shown in Sec. IV that the action of the 
vacuum radiation pressure on the motion of the mirror is a very small perturbation (weak coupling limit). Thus 
the time evolution may be computed from a simple 'semi-classical' model, in which the field evolution is obtained 
assuming the classical prescribed motion 

q{t) = q{0)cos{uJot), (Al) 

where q(t) is the position of the mirror at time t. The dynamical Casimir effect is described by the interaction 
Hamiltonian (see Ref. ||l^, and compare with the first term in Eq. (^)) 

Hint = -m'P- (A2) 
The amplitude b{/St) for the creation of photon pairs corresponding to frequencies uji and 102 at time is given by 

6(o.i,W2;Ai) = J(a;i,W2|7'|0) f^' dt'e'^'^'+^'^'' q{t'). (A3) 
^ Jo 

According to Eq. ( |A3| ) , the amplitude depends on the sign of q, which is very important to the discussion of decoher- 
ence. 



Replacing Eq. (Al) into ([A3|), we find for the two-photon probabilities 



|6(c.i,c.2; Ai)|2 « l^\{0\r\cu,,u;2)M0f (A4) 
n 



^ sin^ [{uji +UJ2- a;o)Ai/2] 



For cjoAi ^ 1, the r.-h.-s. of Eq. ( A4) is sharply peaked around uji + uj2 — wq. Thus, for large times energy is well 



defined, in agreement with the time-energy uncertainty relation. In this limit, Eq. (A4) yields 



|6(c.i,c.2;Ai)|2 « ^mv\u,,0J2)?m^ ^ (A5) 

X(5(ljJi -I- ljJ2 — ^0) 

Since the source of the radiated energy is the motion of mirror, one may expect that the two-photon probabilities are 
related to the amplitude decay rate F. The radiated energy during Ai is 

Y /"OO /"OO 

AE = - duji dw2|6(wi,w2; Ai)|2?i(wi -I-W2), (A6) 
2 Jo Jo 

which according to Eq. ( |A5| ) is proportional to the time interval Z\i. The energy of the mirror decays as dEm/dt = 
— 2F£'Af , where Em — Mg(0)^/2. Hence, from energy conservation we have 

Mq{Oy A£ ' 

leading, with the help of ( [A5| ) and (|A6|), to the representation given by (|4C|). 

To analyze the effect of decoherence, we take the field to be initially in the 'even' superposition state |'0)o = 
(la) -I- I — a))/V^, so that the mirror-plus-field state at i = is 

|*)o - |^)c® |0). 



By linearity, its time evolution is obtained from the two-photon amplitudes given by (A3): 

I*) A* = ® |</'+)At + \-a)® \(j)-)^)/V2, (A7) 

where 
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2^ /'CXD /'OO 

^ Jo JQ 



(A8) 



The already noted sensitivity of the two-photon amphtudes to the phase of the motion of the mirror, which is exphcit 
through the 'minus' sign for \4>^) i n (|A8| ), generates entanglement be twe en mir ror and field. This is discussed in 
Sec. V, whose starting point is Eq. ([f5[), which is derived by replacing (AS) into (A7). 
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